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Abstract 

Motivated by set estimation problems, we consider three closely related shape con- 
ditions for compact sets: positive reach, r-convexity and rolling condition. First, the 
relations between these shape conditions are analyzed. Second, we obtain for the esti- 
mation of sets fulfilling a rolling condition a result of "full consistency" (i.e., consistency 
with respect to the Hausdorff metric for the target set and for its boundary). Third, 
the class of uniformly bounded compact sets whose reach is not smaller than a given 
constant r is shown to be a P-uniformity class (in Billingsley and Tops0e's (1967) 
sense) and, in particular, a Glivcnko-Cantelli class. Fourth, under broad conditions, 
the r-convex hull of the sample is proved to be a fully consistent estimator of an r- 
convex support in the two-dimensional case. Moreover, its boundary length is shown to 
converge (a.s.) to that of the underlying support. Fifth, the above results are applied 
to get new consistency statements for level set estimators based on the excess mass 
methodology (Polonik, 1995). 
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estimation; boundary length; excess mass. 

AMS 2000 subject classifications: Primary 47N30; secondary 60D05, 62G05. 

1 Introduction 

Three geometric closely related properties, with fairly intuitive interpretations, are an- 
alyzed in this paper. They are called positive reach, (see, e.g., Federer [13], Rataj [27] , 
Ambrosio et al. [T] ) , r-convexity, (Perkal [23] , Mani-Levitska [TTj , Walther p2] ) and rolling 
condition (Walther [32l[33]). Our interest in these "rolling- type properties" is motivated 
in the framework of set estimation, that is, the problem of reconstructing a set S C M.'^ 
(typically a density support or a density level set) from a random sample of points; see 
e.g. Cuevas and Fraiman [8] for a recent survey. The classical theory of this subject, 
dating back to the 1960's, is largely concerned with the assumption that S is convex; see 
e.g. Diimbgen and Walther [12j or Reitzner [28] for a survey. The rolling-type properties 
have been employed as shape restrictions on S alternative to (and much broader than) the 
convexity assumption. 
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The rolling-type conditions are useful in statistics and stochastic geometry at least in 
two ways. First, they can be sometimes incorporated to the estimator: for example, if 
the (compact) support S" of a random variable X is assumed to be r-convex one could 
estimate S, from a random sample of X by taking the r-convex hull of the sample points, 
much in the same way as the convex hull has been used to estimate a convex support; 
see Rodriguez-Casal [29] . Second, they can be used as regularity assumptions in order to 
get faster rates of convergence for the estimators; see, e.g., Cuevas and Rodriguez-Casal 
[To] . Pateiro-Lopez and Rodriguez-Casal [22]. Of course other, perhaps more standard, 
regularity conditions have also been used in set estimation. They rely on usual smoothness 
assumptions on the boundary or the underlying density, defined in terms of derivatives. 
See Biau et al. [31 Sj and Mason and Polonik [18] for some recent interesting examples. 
Some deep results on the connection between differentiability assumptions and rolling-type 
conditions can be found in Federer [13] . Walther [33] and Ambrosio et al. pL]. 

The most popular among the rolling-type properties is by far the positive reach con- 
dition. It was introduced by Federer [13] in a celebrated paper which could be considered 
as a landmark in geometric measure theory. Among other relevant results, Federer [13] 
proved that (Theorem 5.6), for a compact set S with reach r > 0, the volume (Lebesgue 
measure) of the e-parallel set B{S,e) can be expressed as a polynomial in e, of degree d, 
for < e < r. This is a partial generalization of the classical Steiner formula that shows 
this property for convex sets for all e > 0. So the positive reach property can be seen as 
a natural generalization of convexity in a much deeper sense than that suggested by the 
definition. For some recent interesting contributions on this property see Ambrosio et al. 
[1] and Colesanti and Manselli [7]. An application in set estimation, more specifically in 
the problem of estimating the boundary measure, can be found in Cuevas et al. [9]. 



The r-convexity property provides a different but closely related generalization of con- 
vexity (see Section [2] for precise definitions). Whereas this property has also a sound 
intuitive motivation it is much less popular. An earlier reference is Perkal [23] but, to our 
knowledge, the first statistical application is due to Walther [32] who uses this condition 
in the setting of level set estimation. A study of the r-convex hull as an estimator of an 
r-convex support can be found in Rodn'guez-Casal [29] . 

Let us know establish some notation and basic definitions. We are concerned here 
with subsets of M'^ although some concepts and results can be stated, with little additional 
effort, in the broader setup of metric spaces. The Euclidean norm in M'* will be denoted 
by II • II . Given a set A C M"^, we denote by A'^, int(^) and dA the complement, interior 
and boundary of A, respectively. We denote by B(x,r) the closed ball with centre x and 
radius r. By convenience, the open ball 'mt{B(x,r)) will be denoted by B{x,r). 

In the problem of estimating a compact set and/or its boundary we need to use some 
suitable distances in order to assess the quality of the estimation and establish asymptotic 
results (concerning consistency, convergence rates and asymptotic distribution). The most 
usual distances in this setting are the HausdorfF distance dn and the distance in naeasure 
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d^. The definitions are as follows. 

Let M be the class of closed, bounded, nonempty subsets of M'^. For A,C C M, the 
Hausdorff distance between A and C is defined by 

dniA, C) = inf {e > : C C B{A, e) and A C B{C, e)} , 

where B{A, e) denotes the closed e-parallel set of A, B{A, e) = {x G M'^ : 6a{x) < e}, with 
5a{x) = inf{|[x — y|| : y £ A}. (M, dj^) is a complete locally compact metric space. 

Let u he a Borel measure on M'^, with i'{C) < oo for any compact C. Let A, C be Borel 
sets with finite z^-measure. The distance in measure between A and C is defined by 

du(A,C) = i/(AAC), 

where A denotes the symmetric difference between A and C, that is, A/\C = {A\C)\J{C\ 
A). Often V is either a probability measure or the Lebesgue measure on R'^, which we will 
denote by /i. Note that the distance function dy is actually a pseudometric but it becomes 
a true metric if we identify two sets differing in a z^-null set. This amounts to work in the 
quotient space associated with the corresponding equivalence relation. 

Besides the estimation of the set S and the boundary 95, some functionals of S are 
sometimes of interest as estimation targets. This is the case of the {d — l)-dimensional 
boundary measure, L{S) (that is the perimeter of for d = 2 and the surface area for 
d = 3). There are several, not always equivalent, definitions for L{S) (see Mattila [H]) but 
we will mainly use the outer Minkowski content, 

L(S) = lim^^i^^MA^. (1) 

Under regularity conditions, the limit in ([T]) coincides with the usual Minkowski content 
given by 

Lo(5)^iim^»». (2) 

e-i-O 2e 

We refer to Ambrosio et al. [1] for general conditions ensuring the existence of the outer 
Minkowski content and its relation to other measurements of the boundary of S. 

The contributions in this paper can be now summarized as follows. In Section [2] we 
clarify the relations between the rolling-type conditions pointing out that if the reach of 
a set is r, then it is r-convex and, in turn, r-convexity entails the r-rolling property. See 
Propositions [J and [21 It is also shown that the converse implications are not true in general. 
In Section [3] we prove that, under very general conditions on and S, (i//(5„,5) — )• 
plus dH{dSn,dS) — )• (we call this simultaneous convergence full convergence) implies 
du{Sn, S") — )■ (Theorem [2]) . We also show (Theorem [3]) that if the sets Sn fulfill the rolling 
condition, the dj^-convergence dniSmS) — ?> implies full convergence. This result (which 
has some independent interest) is used below in the paper. In Section[4]we show that, under 
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broad conditions, the class of uniformly bounded sets with reach > r is a P-uniformity 
class (in the sense of Billingsley and Tops0e [S]) and therefore also a Glivenko-Cantelli class. 
This is interesting from the point of view of empirical processes, see e.g., Devroye et al. 
pT] or van der Vaart [31], and will be also used in Section [5] of this paper. In subsection l5.ll 
we show that, if a support S in is assumed to be r-convex then, the natural estimator 
(which, as mentioned above, is the r-convex hull of the sample) is consistent in all the 
usual senses. In particular, it provides also a plug-in consistent estimator of the boundary 
length whose practical performance is checked through some numerical comparisons in the 
appendix. This is an interesting contribution to the theory of nonparametric boundary 
estimation which so far relies mostly on the use of two samples (one inside and the other 
outside the set S); see Cuevas et al. [9], Pateiro-Lopez and Rodn'guez-Casal [22] and 
Jimenez and Yukich [16] . The results mentioned in the above points are used in subsection 
15.21 for the problem of estimating density level sets of type {/ > A} using the excess mass 
approach (see Polonik [25]). In particular, we obtain new consistency properties for the 
estimation of {/ > A} as well as uniform consistency results (with respect to A) for the 
same problem. 

2 Rolling-type assumptions 

Convexity is a natural geometrical restriction that arises in many fields. While the study of 
this topic dates back to antiquity, most important contributions and applications date from 
the 19th and 20th centuries. We refer to the Handbook edited by Gruber and Wills [14] 
for a complete survey of convex geometry and its relations to other areas of mathematics. 
However, convexity leaves out usual features of sets such as holes or notches and may be 
considered an unrealistic assumption in some situations. As a natural consequence, there 
are meaningful extensions of the notion of convexity; see e.g. Mani-Levitska [17]. In this 
Section, we shall consider three different shape restrictions that generalize that of convexity. 
We group them under the name of rolling-type conditions. Their formal definitions are as 
follows. 

Definition 1. Given r > 0, a set S C M.'^ is said to fulfill the (outside) r-rolling condition if 
for all X € dS there is a closed ball with radius r, Bx, such that x S Bx and int{Bx)r\S = 0. 

The radius r of the ball acts a smoothing parameter. We refer to Walther [33] for a 
deep study of this condition. 

Following the notation by Federer [13], let Unp(S') be the set of points x having 
a unique projection on S, denoted by is{x)- That is, for x G Unp(S'), ^^(x) is the unique 
point which minimizes 6s{x). 

Definition 2. For x ^ S let reach(S', x) = sup{r > : B{x,r) C Unp(5')}. The reach of 
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S is then defined by 

reach(5) = mf{reach(S', x) : x ^ S} 
and S is said to be of positive reach i/reach(S') > 0. 

Besides all convex closed sets (which have infinite reach) and regular submanifolds of 
class 2 of M'^, the class of sets of positive reach also contains nonconvex sets or sets whose 
boundary is not a smooth manifold. Sets with positive reach were introduced by Federer 
[13] who also obtained their main properties. In particular, these sets obey a Steiner 
formula in the following sense. 

Theorem 1. (Federer (1959, Th. 5.6)) Let S C M."^ a compact set with reach(S') > 
0. Let K be a Borel subset of . Then there exist unique Radon measures '^q{S, ■), 
■),..., ^diS, •) over R'^ such that, /or < e < reach(S), 

d 

fi{B{S, e)n{x: ^s{x) G K}) = J] e^-'bd^i^S, K) (3) 

1=0 

where bo = \ and, for j > 1, bj is the j-dimensional measure of a unit ball in . 

Remark 1. Theorem {1\ plays an important role in some problems of stochastic geometry 
regarding the calculation of boundary measurements. The measures := ^i{S,-) are 
the so-called curvature measures associated with S, for j = 0,1, . . . ,d. In particular, if 
reach(5) > 0, it is straightforward from Theorem [T] with K = S, that the outer Minkowski 
content L{S) is finite, and corresponds to the first order coefficient in the expansion 
Furthermore, by Corollary 3 and inequality (27) in Ambrosio et al. [I], combined with 
Federer's Theorem, we have that if reach(S') > 0, then dS has finite Minkowski content 
Lq{S), which in particular yields nidS) = 0. For this reason we do not impose in what 
follows fj,{dS) = whenever we assume positive reach for S. 

Definition 3. A set S C is said to be r- convex for some r > if S = Cr{S), where 

B{y,r)n5=0 

That is, S is r-convex if for all x € S^, there exists y € B{x, r) such that B{y, r) n S" = 
0. We refer to Perkal [23] for elementary properties of r-convex sets and connections 
between convexity and r-convexity. Walther [33] and Rodrfguez-Casal |_29J also deal with 
this shape restriction in the context of set estimation. A natural question is whether certain 
characterizations of convex sets are still meaningful in the context of generalized notions of 
convexity. For instance, given a set S and r > it is possible to find the minimal r-convex 
set containing S, the so-called r-convex hull of S. In fact, it follows from the properties 
of r-convex sets that the r-convex hull of S coincides with the set Cr{S) given in ([4]), see 
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Perkal [23]. Note that the definition of Cr{S) resembles that of the convex huh (with bahs 
of radius r instead of halfspaces). However, the same property does not hold when we 
consider the reach condition. It is not always possible to define the so-called r-hull of S, 
that is, the minimal set containing S and having reach > r, see Colesanti and Manselli [7]. 
In fact, when S admits such a minimal set it happens to coincide with Cr{S), see Corollary 
4.7 by Colesanti and Manselli f^. This result provides an indirect proof of Proposition [T] 
below, that states that every set with reach r is also r-convex. Proposition [2] establishes 
the relation between r-convexity and the rolling condition. 

Proposition 1. Let S CM.'^ be a compact set with reach(S') > r > 0. Then S is r-convex. 

Proof. Let x £ 5^. If Ss{x) > r, then B{x, r) n 5 = 0. If Ss{x) < r, let S,s{x) be the unique 
point in S such that 6s{x) = \\x — ^s{x)\\. Then 

where Nor(5', (,s{x)) is the set of all normal vectors of S at ^s(x). Define y\ = ^^(x) + \r]x 
with < A < r and take A € (0,r). Now, reach(5, ^5(rc)) > r > A and by part (12) of 
Theorem 4.8 in Federer [13], we get that 7]^ = X~^v where = A and £,s{Csix) + v) = 
^s{x). That is, Csiux) = Csi^)- It is clear that x G B{yr,r). Moreover, B{yr,r) n = 0. 
Indeed, if z G B{yr,r), then z € B{yx,X) for some yx = Csix) + Xrjx with < A < r and 
B{yx, A) n 5 = as a consequence of Csiux) = Cs{x) and reach(S', ^s'(x)) > r > A. □ 

Remark 2. Borsuk's conjecture on local contractibility of r-convex sets. The converse of 
Proposition [1] is not true in general, see Figured] (a) for an example of a r-convex set that 
has not reach r. Even so, we prove in Theorem [6] that if S is a compact r-convex support 
in fulfilling a mild regularity condition (which we call interior local connectivity; see 
Section [5] for details) then S has positive reach, though not necessarily r. If we do not 
assume any additional regularity condition on S, then the conclusion of positive reach does 
not seem that simple to get. This is closely related to an unsolved conjecture by K. Borsuk 
(see Perkal [23] and Mani-Levitska |17j): Is an r-convex set locally contractihle? Note that 
proving that a compact r-convex set has positive reach would give a positive answer to 
Borsuk's conjecture since, according to Remark 4.15 in Federer [13], any set with positive 
reach is locally contractible. Recall that a topological space is said to be contractihle if it 
is homotopy equivalent to one point. In intuitive terms this means that the space can be 
continuously shrunk to one point. The space is called locally contractible if every point 
has a local base of contractible neighborhoods. 

Proposition 2. Let S C M*^ be a compact r-convex set for some r > 0. Then S fulfills the 
r-roUing condition. 

Proof. Let x € dS. Since x is a limit point of the set S'^, there is a sequence of points 
{xn}, where x„ G S'^, that converges to x. By the r-convexity, for each n there exists 
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Hn G B(xn,r) such that B{yn,r) fl 5 = 0. Since {?/«} is bounded it contains a convergent 
subsequence which we denote by {un} again. Then y-n ^ y and it is not difficult to prove 
that B{y,r) H 5 = and — < r. Now, a; € 5 since S is closed and, therefore, 
X € dB{y,r), which concludes the proof. □ 

The converse implication is not true in general, see Figure [T](b) for an example of a set 
fulfilling the r-rolling condition but not r-convex. 



(a) (b) 




Figure 1: (a) r-convex ^ r-reach. (b) r-rolling ^ r-convex. 



3 Boundary convergence and full convergence in sequences 
of sets 

As mentioned in the introduction, the focus in this paper is on the reconstruction (in 
the statistical sense) of an unknown support S from a sequence of estimators {5„} based 
on sample information. In many practical instances, including image analysis, the most 
important aspect of the target set is the boundary dS. However, it is clear that, even in 
very simple cases, the Hausdorff convergence dniSmS) — )• does not entail the boundary 
convergence dH{dSn,dS) 0. See e.g. Bafllo and Cuevas [2\, Cuevas and Rodriguez-Casal 
[To] and Rodn'guez-Casal [29] for some results on boundary estimation. We introduce the 
following notion of convergence. 

Definition 4. Let {Sn} be a sequence of compact non-empty sets in W^. Let S C be a 

compact non-empty set. We say that we have full convergence of{Sn} to S ifdniSn, 5) — )• 
and dnidSn, OS) 0. 

The following result shows that, under very general conditions, the Hausdorff conver- 
gence of the sets and their boundaries implies also the convergence with respect to the 
distance in measure. This accounts for the term "full convergence". 
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Theorem 2. Let {Sn} he a sequence of compact non-empty sets in endowed with a 
Borel measure v , with ^{C) < oo for any compact C . Let S he a compact non-empty set 
such that v{dS) = 0, dnlSn, S) ^ and dnidSn, dS) 0. Then d^{Sn, S) 0. 

Proof Take e > 0. We first prove that, for n large enough, 

{xeS: 59s{x) > 2e} C Sn (5) 

and 

{xeSn: SqsAx) > 2e} C S. (6) 

To see ([5]) take x G S such that 5qs{x) > 2e and n large enough such that S C B{Sn-,() 
and dnidSn, dS) < e. If x ^ Sn, then x £ B{Sn, e) \ Sn- Therefore 5qs„{x) < e so that 

6ds{x) = dH{{x},dS) < dH{{x},dSn) + dH{dS,dSn) < 2e, 

which contradicts 5qs{x) > 2e. 

The proof of ([6]) follows along the same lines. 
Finally, note that as a consequence of ([S]) and 

SnAS c{xgS: 69s{x) < 2e} U {x G 5„ : SasAx) < 2e}, 

which for large enough n is a subset of {x G S : Sgsix) < 3e}, that decreases to dS as 
e I 0. Since u is finite on bounded sets and u{dS) = this entails limsup z^(S'„A5) = 
which concludes the proof. □ 

We next show that the above result applies to the important class of sets fulfilling the 
r-rolling condition. Therefore, as shown in Propositions [1] and [21 it also applies to the class 
of r-convex sets and to that of sets with positive reach. 

Theorem 3. Let {Sn} be a sequence of compact non-empty sets in satisfying the r- 
rolling condition. Let S he a compact non-empty set such that dniSmS) — > 0. Then, 
(a) dH{dSn,dS) ^0. 

(h) //reach(S'„) > r > 0, then we also have du{Sn-,S) — )• for any Borel measure v 
(finite on compacts) absolutely continuous with respect to the Lehesgue measure /U. 

Proof, (a) Assume that the result is not true. Then, 

(i) There exists e > such that for infinitely many n € N there exists x„ G dS with 

^dsA^n) > e or 

(ii) There exists e > such that for infinitely many n G N there exists Xn G dSn with 
^ds{xn) > e- 
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First, assume that (i) is satisfied. Since 5 is compact, there exists a convergent sub- 
sequence of {xn} which we will denote again as the original sequence. Let x G dS be the 
limit of {xn}- The Hausdorff convergence of Sn to S implies that Ss„{x) < e/2 for infinitely 
many n. Furthermore, for large enough n, 

^dSni^) > ^dsA^n) - \\x -Xn\\ > e/2. 

Now, 5s„{x) < e/2 together with 5qs„{x) > e/2 yields x G int(S'„) for infinitely many n 
and B{x,e/2) n S*^ = 0. But, since x G dS and S is closed, we can consider y & such 
that ||x — y\\ < e/2 and 5s{y) > 0. Again by the Hausdorff convergence of Sn to S we get 
that y (z for infinitely many n which yields a contradiction. 

Assume now that (ii) is satisfied. We can assume e < r. The Hausdorff convergence of 
Sn to S implies that 5s{xn) < e for infinitely many n. This, together with 5Qs{xn) > e, 
yields x„ G int(S') for infinitely many n. Let < A < e/2 < r. Since Xn G dSn and 
the sets 5„ satisfy the r-rolling property, there exists for each n G N a ball B{cn,\) such 
that Xn G dB{cn-,X) and B{cn-,X) PI 5™ = 0. Again, let us denote by {xn} a convergent 
subsequence of {x^}. Then with X G S and 5qs{x) > e/2, that is, B{x,e/2) C S. 

Now, let {cn} be a convergent subsequence of {c„}. Then c„ — t- c, with c G S"^ and 
Ssr,ic) > A/2 for infinitely many n. By the Hausdorff convergence of Sn to 5, we get that 
c & S^ which yields a contradiction with B{x, e/2) C S. 

(b) By Theorem 4.13 in Federer [I3], reach(5) > r and, in particular, fJ-{dS) = 0, see 
Remark [TJ The result is now straightforward from (a) and Theorem [2j □ 

Remark 3. If A is a class of sets, closed with respect to the Hausdorff topology, and the 
dn -convergence in this class implies the dn -convergence of the respective boundaries, then 
Theoreni\^ entails that this class is also closed with respect to d^ whenever u^dA) = for 
all A ^ A. So, the du -compactness of A implies the dy- compactness of this class. 

On the other hand. Theorem [3 shows that this in fact applies to the class of sets with 
reach > r, which is dn -closed from Theorem 4-13 and Remark 4-14 Federer ]13^ . Thus, 
the class of subsets with reach > r of a compact set is d^-compact. This will be useful below 
(subsection \5.S^) in order to apply the results in Polonik f25\j for classes of sets defined in 
terms of reach properties. 

4 P-uniformity: Billingsley-Tops0e theory and its applica- 
tion to classes of sets with positive reach 

Let Xi , X2 , . . . be a sequence of independent and identically distributed random elements 
on a probability space (O, 3", P) with values in a measurable space {E, 23), where E is 
a metric space and S stands for the Borel cr-algebra on E. Denote by P the probability 
distribution of Xi on S and let P„, be the empirical probability measure associated with 
^1 ) • • • ) Xn ■ 
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The almost sure pointwise convergence on H of P„ to P is ensured by the strong law 
of large numbers. Moreover, for appropriate classes A C "B the uniform convergence 

sup \Fn{A) - P{A)\ 0, a.s. (7) 

also holds. 

A class A of sets fulfilling ([7]) is called a Glivenko-Cantelli class (GC-class). They are 
named after the classical Glivenko-Cantelli theorem which establishes the result ([7|) for the 
case where E = U. and A is the class of closed half lines A = (— oo,x], x G M. The study 
of uniform results of type ([7]) is a classical topic in statistics. A well-known reference is 
Pollard pi]. A useful summary, targeted to the most usual applications in statistics, can 
be found in Chapter 19 of van der Vaart |31j . 

A popular methodology to obtain Glivenko-Cantelli classes is based on the use of the 
well-known Vapnik-Cervonenkis inequality (see, e.g. Devroye et al. [H]) which relies on 
combinatorial tools. This inequality provides upper bounds for lP{sup^gyi |Pn(^) — P(A)| > 
e} which depend on the so-called shatter coefRcients and VC-dimension of the class A. 
However, this approach is not useful in those cases where the VC-dimension of 71 is infinite. 
This is the case where A is the class of closed convex sets in and therefore also for the 
class of closed r-convex sets. 

Nevertheless, it can be shown that the family of all convex sets in is a GC-class. This 
can be done following an alternative approach (maybe less popular than the VC-method), 
due to Billingsley and Tops0e [5] ; see also Bickel and Millar pj . This methodology is rather 
based on geometrical and topological ideas and therefore turns out to be more suitable for 
set estimation purposes. The basic ideas of Billingsley- Tops0e approach can be summarized 
as follows. 

Let Pn and P be probability measures on 25. A set A is called a P-continuity set if 
P{dA) = 0. The sequence Pn is said to converge weakly to P if 

Pn{A)^P{A) 

for each P-continuity set A C "B. A subclass A C H is said to be a P-continuity class if 
every set in 71 is a P-continuity set. A subclass 71 C 23 is said to be a P-uniformity class if 

sup\Pn{A)-P{A)\^0 
AeA 

holds for every sequence Pn that converges weakly to P. Note that the P-uniformity 
concept is not established just for sequences of empirical distributions but in general for 
any sequence of probability measures converging to P. Billingsley and Tops0e [5j derived 
several results establishing conditions for a class 71 to be a P-uniformity class. As a 
consequence, some useful criteria for obtaining GC-classes immediately follow. 

The following theorem provides three sufficient conditions to ensure that a class 71 is a 
P-uniformity class. 
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Theorem 4. (Billingsley and Tops0e (1967, Th. 4)) If E is locally connected and if A is a 
P-continuity class of subsets of E, then each of the following three conditions is sufficient 
for A to he a P -uniformity class. 

(i) The class dA = {dA : A € ^1} is a compact subset of the space M of non-empty 
closed bounded subsets of E. 

(a) There exists a sequence {C„} of bounded sets with P(int(C„)) — t- 1 and such that, for 
each n, the class d{Cn r\A) = {d{Cri r\ A) : AG A} is a compact subset of M. 

(Hi) There exists a sequence {C„,} of closed, bounded sets with P(int(C„)) — t- 1 and such 
that, for each n, the class Cn H dA = {C„ fl dA : A € A} is a compact subset of M. 

Theorem S] can be used to prove that the class of all convex sets in is a P-uniformity 
class, and in particular, a Glivenko-Cantelli class. The following theorem provides a partial 
extension of this property: we show that, under an additional boundedness assumption, 
convexity can be replaced with the broader condition of having a given positive reach. 
Again, the basic tool is Theorem HI 

Theorem 5. Let K be a compact non-empty subset of and A = {A C K : A ^ 
0, A is closed, and reach(A) > r}, for r > 0. Then the class dA = {dA : A € A} is a 
compact subset o/M. Moreover, A is a P-uniformity class for every probability measure 
P such that P is absolutely continuous with respect to the Lebesgue measure /i. 

Proof. Let {dAn} be a convergent sequence of sets in dA. By Federer's closeness theorem 
for sets of positive reach (Theorem 4.13 in Federer [13]) it follows that the class A is compact 
with respect to the Hausdorff metric and, therefore. An has a convergent subsequence whose 
limit is a set A € A. In view of Propositions [T] and [21 we can apply Theorem [3] to the 
class A, and this yields dAn — ?• dA in dn. Now, from Theorem [J] (i), A is a P-uniformity 
class. □ 

5 Applications 

5.1 Estimation of r convex sets and their boundary lengths 

Estimation of r- convex supports 

As indicated in Section [21 r-convexity is a natural extension of the notion of convexity. 
Prom the point of view of set estimation, r-convexity is particularly attractive as the 
estimation of an r-convex compact support S from a random sample Xi, . . . ,Xn drawn 
on S, can be handled very much in the same way as the case where S is convex. In 
this classical situation (which has been extensively considered in the literature, see, e.g., 
Diimbgen and Walther J12]) the natural estimator of S is the convex hull of the sample. 
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In an analogous way, if 5 is assumed to be r-convex, the obvious estimator is the r-convex 
huh of the sample, which we will denote by 5„. Recall from Q that 

Sn= n B{y,rr. (8) 

S(y,r)n{Xi,...,X„}=0 

This estimator can be explicitly calculated in a computationally efficient way (at least 
in the two-dimensional case) through the R-package alphahull; see Pateiro-Lopez and 
Rodriguez-Casal [2T]. In Figure [2] we show an example of the r-convex hull estimator. 
Note that the boundary of the r-convex hull estimator is formed by arcs of balls of radius 
r (besides possible isolated sample points). The arcs are determined by the intersections 
of some of the empty balls that define the complement of the r-convex hull, see Equation 




Figure 2: In gray, r-convex hull of a sample of size n = 200 from a uniform distribution on 
the set S = B{0, 0.5) \ B{0, 0.25). The value of r is 0.25. 



The r-convex hull estimator was first considered, from a statistical point of view, in 
Walther [32]. Rodn'guez-Casal [29] studied its properties as an estimator of S and dS, 
providing (under some rolling- type assumptions for S) convergence rates for dniSmS), 
dH{dSn,dS) and d^{Sn,S) which essentially coincide with those given by Diimbgen and 
Walther p!2] for the convex case. 

Estimation of the boundary length 

We are interested in estimating the boundary length of S, L{S), as defined by the outer 
Minkowski content given in ([1]). Recall that (see Remark [1]) if reach(S') > 0, then the outer 
Minkowski content is finite. 
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In Pateiro-Lopez and Rodriguez- Casal [22] the estimation of the usual (two-sided) 
Minkowski content Lq{S) given in ^ is considered under double smoothness assump- 
tions of rolling-type on S (from inside and outside dS) which in fact are stronger than 
r-convexity. Under these assumptions, these authors improve the convergence rates for the 
Minkowski content Lq(S) obtained in Cuevas et al. [9]. Another recent contribution to 
the problem of estimating boundary measures is due to Jimenez and Yukich [16]. In all 
these cases the estimators are based on a sample model which requires random observations 
inside and outside S in such a way that for each observation one is able to decide (with no 
error) whether or not it belongs to S. 

Theorem [6] provides further insights on the approach of these papers in the sense that, 
under minimal assumptions, gives a fully consistent estimator Sn of an r-convex set 5 C 
and a plug-in consistent estimator L{Sn) of L{S) based on a unique inside sample. The 
performance of L{Sn), in terms of bias and variance is analyzed in the appendix. 

Before stating Theorem [6] we need some preliminaries. 

Definition 5. We will say that a set S fulfills the property of interior local connectivity 
(ILC) if there exists ao > such that for all a < uq and for all x ^ S, int(i?(x, a) fl S) is 
a non-empty connected set. 




Figure 3: The set S in gray does not fulfill the ILC property at the central point x. 



Let us define an r -circular triangle as a compact plane figure limited by three sides: 
two of them are arcs of intersecting circumferences of radius r; the third side is a linear 
segment which is tangent to both circumference arcs. In other words, an r-circular triangle 
is an isosceles triangle with a linear segment as a basis and two r-circumference segments 
as the other sides. See Figure [H 

Theorem 6. Let S C M? be a compact r-convex set with fJ-{S) > that fulhlls the ILC 
property. Let Xi, . . . ,Xn he a random sample drawn from a uniform distribution with 
support S. Denote by Sn the r-convex hull of this sample, as defined in Then, 

(a) Sn is a fully consistent estimator of S. 

(b) S has positive reach and L{Sn) — ^ L{S), a.s., L{S) being the outer Minkowski 
content of S. 



13 




Figure 4: In gray, an r-circular triangle. 



Proof, (a) Since {Xi,...,Xn} C 5„ C S and dndXi, . . . , Xn}, S) — )• 0, a.s. we also 
have dniSnjS) — )• 0, a.s. In order to get the full consistency we only need to prove 
dnidSmdS) 0, a.s. but this follows directly from Proposition [2] and Theorem [3l 

(b) Define Sn = Sn \ I{Sn), where 1(5^) is the set of isolated points of 5„. Thus, for 
large enough n, with probability one Sn is a compact (not necessarily connected) set whose 
boundary is the union of a finite sequence of r-arcs, that is, circumference arcs of radius 
r. We will say that a boundary point of Sn is "extreme" if it is the intersection of two 
boundary r-arcs from different r-circumferences. The proof is based on the following two 
lemmas. 

Lemma 1. Under the conditions of Theorem\^ there exists tq > such that reach(S'„) > rg 
a.s. for all n. 

Proof. If the sequence reach(S'„) were not bounded from below a.s. we would have, with 
positive probability, a sequence of points with z„ ^ Sn and another sequence {xn,yn), 
where x„ and y„ are boundary points of Sn, with x„ ^ y„ and such that 

\\Zn - Xn\\ = \\Zn - UnW = mf \\Zn - x\\ = r„ 0. (9) 

In principle, the "projection points" Xn and yn could be extreme points (as defined above) 
or "boundary inside points", that is, points belonging to a boundary r-arc but different 
from the extremes of this arc. However, it is easily seen (Figure[5]) that the latter possibility 
(i.e., Xn or y„ is a boundary inside point) leads to a contradiction: Indeed, note that if 
Xn belongs to an r-arc in dSn with extremes a and b and Xn a, Xn ^ b, then at a fixed 
distance rn < r there is only one point in the complement of Sn (necessarily equal to Zn) 
such that Xn is the projection of Zn on Sn- This would imply that the projection of z„ is 
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Figure 5: If Zn has two projections x„, and ?/„ onto 5„, neither x„ nor y„ can be boundary 
inside points. 



unique since the arc with extremes a and b is an arc of a circle of radius r that does not 
intersect 5„. 

Thus, Xn and yn must be extreme points and our proof reduces to see that we cannot 
have do]) with a positive probabihty. So, let us assume that ([9]) were true with a positive 
probability (let us denote by A the corresponding event). 

Let and be two circular triangles, with vertices x„ and respectively, deter- 
mined by the r-arcs in dSn whose intersections are Xn and yn, see Figure [H Note that 
these triangles are not necessarily included in Sn] just a portion of each triangle close to 
the vertex is in general included in Sn- 

First, let us prove that the "heights" (i.e. the distances from the vertices to the bases) 
of and must necessarily be bounded from below a.s. on A. To see this, recall that 
each Xn is the intersection of the boundary arcs of two balls of radius r, B(cin,r) and 
B{c2n,r), whose interiors do not intersect Sn- By construction, height(r^) — >■ for some 
subsequence implies that the centers (cin,C2n) must fulfill ci = limci^ = limc2n = C2. 
This means that the r-circumferences providing the boundary of Sn at both sides of Xn 
tend to coincide as n tends to infinity. We will prove that ci = C2 leads to a contradiction. 
Indeed, if ci = C2 we would have that, for large enough n, the centers ci„ and C2n would be 
very close and the boundary arcs, dB{cin,r) and dB{c2n,r), would intersect to each other 
not only at x„ but also at another point x* such that ||x„ — x* || > 7 > 0. We conclude 
B{zn,rn) \ {xn} C int(i?(ci„, r) U B{c2n,r)) C Sn , for large enough n (to see this note 
that the boundary of B{cin,r) U B{c2n-,i') near x„ coincides with the triangular sides of 
T^). Now, from ([9]) we have obtained a contradiction ?/„ G dB{zn,rn) n Sn and y„ 7^ x„. 

As the space of compact sets endowed with the Hausdorff metric is locally compact, 
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Figure 6: In gray, 



Tn and T^. 



there exist a.s. convergent subsequences of and Z!^, {xn} and {y^} which we will denote 
again as the original sequences. 

Let and be the a.s. Hausdorff limits of and T^, respectively. Since the heights 
of and are bounded from below, and must be also non-degenerate circular 
triangles. Denote x = lima;„ = limy„ = limz„, a.s. Note that the fact that z„ ^ 5„ 
ensures that x € dS. Then, we have two possibilities, 

(i) if n = {x} we would get a contradiction: to see this note that T = U is 
just the union of two r-circular triangles which are disjoint except for the common vertex 
{x}. Thus, for e small enough B{x, e) fl 5 is included in T. This contradicts the ILC of S 
at X. See Figure [7] (a) . 

(ii) On the other hand, the possibility n 7^ {x} leads also to contradiction. As 
the vertices of and tend to the same point, we would necessarily have that, with 
probability one, yn (or x„) belongs to the interior of one of the r-circles defining the arcs 
of (or T^), see Figure [7] (b). However, by construction, the r-circles defining the arcs 
of and do not intersect Sn ■ This concludes the proof of Lemma [H □ 

Lemma 2. Under the conditions of Theorem\^ 



Proof. Assume that (llOp is not true. Then, with positive probability we have either 

(i) There exists e > such that for infinitely many n € N there exists Xn € Sn with 
Ssixn) > e or 

(ii) There exists e > such that for infinitely many n G N there exists Xn (z S with 



dH{Sn,S) 0, a.s. 



(10) 
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(a) (b) 




Figure 7: (a) Ti nT2 = {x} yields a contradiction with the ILC property, (b) Ti nT2 ^ {x} 
yields a contradiction with the fact that the circles of radius r defining the arcs of and 
cannot intersect Sn- 

First, since S is r-convex, Sn C S with probability one. Therefore Sn C Sn C B{S,e) 
for all e > and we cannot have (i). 

With regard to (ii), let {xn} be a sequence of points in S. Since S is compact, there 
exists a convergent subsequence which we will denote again {xn}- Let x € S he the limit 
of {xn} and y £ int(5) such that — y|| < e/2. Note that the existence of such y follows 
from the ILC property. Then, there exists eo > (eo <^ r) such that -B(y, eo) C S. Using 
(a) in Theorem [6] and a similar reasoning to that of ([5]) we get that, for n large enough, 
B{y,ei) C Sn for some ei < eq. Finally, for n large enough, 

<55„(2;„) < SB(y,ei)iXn) < SB{y,ei)ix) + \\Xn - x\\ < \\x - y\\ + ||x„ - x|| < 6. 

That is, (ii) neither can be true. This concludes the proof of Lemma [21 □ 

Now, in view of Lemmas [1] and [21 the assumptions of Theorem 5.9 in Federer [13] 
are fulfilled (see also Remark 4.14 in that paper). Theorem 5.9 in Federer [13] essen- 
tially establishes that S has positive reach and the curvature measures are continuous 
with respect to dn (see Remark 5.10 in Federer [13j). In particular we obtain that 
^d-iiSn, K) — )• ^d-iiS, K) a.s. for any closed ball K such that S C K. Using Remark 5.8 
in Federer [13] and Sn C S we get that ^d-i{Sn,K) = ^d-iiSn, K ndSn) = ^d-iiSn,dSn) 
and also ^d-i{S, K) = ^d-i{S,dS). The proof of Theorem [6| (b) concludes noting (see 
Remark[II that L{Sn) = L{Sn) = <^d-i{Sn^dSn) and L{S) = <^d-i{S,dS). 

□ 

Remark 4. This result provides, using stochastic methods, a partial converse of Proposi- 
tionUl' we prove that r -convexity implies positive reach for ILC sets in M?. Thus we also get 
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a partial answer to Borsuk's question (is an r -convex set locally contractihle? ) mentioned 
in Remark\^ since from Federer 113], Remark 4-15, any set with positive reach is locally 
contractihle. 

5.2 Applications to the excess mass approach 

Typically, the results on uniform convergence, of Glivenko-Cantelli type, are useful in 
order to establish the consistency of set estimators which are defined as maximizers of 
appropriate functionals of the empirical process. An interesting example in set estimation 
is the excess mass approach, proposed by Hartigan [15j and Miiller and Sawitzki [20] and 
further developed by Polonik [25j and Polonik and Wang [26j, among others. 

The basic ideas of this method can be simply described as follows: given A > 0, denote 
by P an absolutely continuous distribution with respect to the Lebesgue measure ^. Let / 
denote the corresponding /^-density. Define the excess mass functional on the class 23 (M'^) 
of Borel sets, 

Hx{A) = P{A) - Xfi{A), for A G B(R'^). 

Since Hx{A) = f^{f — X)d^, it is obvious that H\{A) is maximized by the level set 
^ = {/> A}. 

This suggests a method to define an estimator of S{\) = {/ > A} which can incorporate 
some shape restrictions previously imposed on this set. Let us assume that ^(A) belongs 
to some given class of sets A (for example, the class of compact convex sets or the class of 
compact r-convex sets in M^). As ^(A) is the maximizer on A of the unknown functional 
Hx{A), we could define ^^(A) as the maximizer on A of the empirical excess mass functional 

Hn,x{A) = ^n{A) - \^i{A), for A A. 

Proposition 3. Let A he the class of compact sets A with reach(A) > r included in a 
given ball B{0,R). Given a sample Xi, . . . ,Xn from an absolutely continuous distribution 
P with density f in W^, let 5„(A) denote the empirical level set estimator defined by 5'„(A) = 
aigmmj^^jiHn.\{A) . Assume that the level set S{X) = {/ > A} belongs to A. Then, 
(a) The estimator Sn{X) is fully consistent to S{X). 

(h) L{Sn{X)) — > L{S{\)), a.s., where L{A) denotes the outer Minkowski content of A. 

Proof, (a) For simplicity, denote i^n,A = -f^m H\ = H, Sn{X) = Sn and ^(A) = S. First, 
note that, by definition, H is a d^-continuous functional and, from Theorem O it is also 
continuous with respect to dn- We have, with probability one, 

I sup HniA) - sup H{A)\ < sup |H„(A) - HiA)\ < sup |P„(A) - P{A)\ ^ (11) 

A£A AeA AgA A£A 

since, from Theorem [5l A is a P-uniformity class. 

From Theorem [3] and Proposition [2] to prove the full consistency we only need to 
establish the d^-convergence 

dHiSn,S)^0, a.s. (12) 
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Now, let us take a value w E 17 ($7 is the common probability space in which the random 
variables Xi are defined) for which du {Sn{^)-, S) — )• does not hold but ([TT]) holds. Then, 
as A is compact, we should have dn {Sn{uj),T) for some subsequence {Sn{uj)} of 
Sn = Sn{^) and for some T ^ S , T & A. Then, (jlip and the continuity of H implies that 
H{T) = H{S). However, the convergence to T 7^ 5 is not possible, since 5 G yi is the 
(unique) maximum of H in A. Thus we should have dn {Sn{^), 5") — ?• for all uj such that 
(|lip holds. We conclude that ()12p must hold with probability one. 

(b) Follows again directly as a consequence of (a) together with Remarks 4.14 and 5.8 
and Theorem 5.9 in Federer [13]. □ 



Our main point in this subsection is to show that our P-uniformity results fit in the 
framework developed by Polonik [25]. This author obtains results of consistency (uniform in 
A) of level set estimators, of type sup;^ (i^(5„(A), ^(A)) — > 0, a.s. One of his key assumptions 
is that the involved classes S are Glivenko-Cantelli classes. To be more specific, the main 
result in Polonik's paper (which is Theorem 3.2) can be applied to our new P-uniformity 
class of sets with reach > r. 

Let us first briefiy recall the basic assumptions for the result in Polonik [25]. 

(Al) For all A > the excess mass functional H{C) = fdfj, — Xfi{C) and its empirical 
counterpart Hn{C) = Pn(C') — ^f^iC) attain their maximum values on the class C at 
some sets denoted by r(A) and r„(A), respectively. 

(A2) The underlying density / is bounded in M'^. 

(A3) e is a GC-class of closed sets with 0^6. 

Note that the maximizer over C, r(A), of H{C) will only coincide with the level set 
S'(A) = {/ > A} whenever 5(A) € C. For this reason r(A) is called the generalized X-cluster. 

Now, let us denote by g* the "measurable cover" of any function g (which of course 
coincides with g when it is measurable). The main result in Polonik's paper is as follows: 

Theorem 7. (Polonik (1995, Th. 3.2)) Let A C [0, 00). Suppose that, in addition to 
(A1)-(A3), the following two conditions hold: 

(i) For a distribution Q in with strictly positive density, the space {Q^dq) is compact. 

(a) For every A € A the generalized X-cluster T(X) is unique up to P-null sets. 

Then, with probability one, 

supdp(r(A),r„(A))* ^0. (13) 

AeA 

Then, our point is that this result applies directly to the case in which C is the class 
A considered in Theorem [5l plus the empty set. This is made explicit in the following 
statement. 
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Theorem 8. Under assumptions(Al), (A2) and (ii) of Theorem^ the conclusion / fi3|) is 
valid if we take Q = AiJ {$}, where A is the class of compact subsets of a given ball B(0, R) 
with reach > r. 

Proof. We only need to prove that conditions (i) and (A3) hold in this case. The validity 
of (i) follows easily since A is d/^-compact from Theorems [3] and O and therefore, from 
Theorem [2l it is also dg-compact for any absolutely continuous distribution Q. Also 
Q = AU {0} is dg-compact. 

As for the condition (A3) it also holds as a direct consequence of Theorem [3 □ 

Note that none of the typical examples of GC-classes (convex sets, ellipsoids, balls,...) 
allows us to consider multimodal densities. In this sense, Theorem [8] can be seen as an 
extension of uniformity results in the excess mass approach beyond the realm of convex 
sets. The class of sets with reach bounded from below is much larger and includes non- 
connected members which are now candidates to be considered as possible level sets for 
multimodal densities. 



Appendix: some numerical comparisons 

In Theorem [6] we have proved that the boundary length L{S) of an r-convex set 5 C 
can be consistently estimated in a plug-in way by L{Sn), where Sn is the r-convex hull of 
the sample. 

Another estimator of L[S) (in fact of Lq{S)) has been recently proposed by Jimenez 
and Yukich [16], based on the use of Delaunay triangulations. This estimator does not rely 
on any assumption of r-convexity but requires the use of sample data inside and outside 
the set S. The numerical results reported by Jimenez and Yukich [I6] show a remarkable 
performance of their sewing based estimator (denoted by L^{S)) which in fact outperforms 
that proposed in Cuevas et al. [9]. 

Our plug-in estimator is not directly comparable with Lfj^{S) since the required sampling 
models are different in both cases. Moreover, L{Sn) incorporates the shape assumption 
of r-convexity on the target set S. Nevertheless, it is still interesting analyzing to what 
extent the use of the r-convexity assumption in L{Sn) could improve the efficiency in the 
estimation. 

We have checked this, through a small simulation study. We have considered two r- 
convex sets S, defined as the domains inside two well-known closed curves: the Catalan's 
trisectrix (as in Jimenez and Yukich [16]) and the astroid. See Figure [H The "true" 
maximal value of r in the first set is r = oo (since it is convex). For the second set r is 
close to 1. In practice, these values are not known so one must assume them as a model 
hypothesis keeping in mind that small values of r correspond to more conservative (safer) 
choices. Recall that the class of r-convex sets is increasing as r decreases. 
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Figure 8: The supports in the simulation study are the domains inside a Catalan's trisectrix 
(left) and an astroid (right). 



The considered sample sizes are n = 5000, 10000. In the case of the estimator in Jimenez 
and Yukich [IB], the sample observations are uniformly generated on a square containing 
the domain S and it is assumed that we know (without error) whether a sample point 
belongs to S or not. In the case of our estimator L{Sn) all the observations are drawn from 
a uniform distribution on S. The results are summarized in Tabled) The reported values 
in the columns "mean" and "std" correspond to averages over 1000 runs. The true values 
of L(S) are 20.7846 (for the Catalan's trisectrix) and 6 (for the astroid). The outputs show 
a better behavior of L(5„), especially in terms of variability. 

Table 1: Mean and standard deviation of the sewing-based estimator (5") and the plug- 
in estimator L{Sn), where Sn is the r-convex hull of the sample. The estimator Sn is 
computed for different values of r. The reported values correspond to averages over 1000 
runs of samples size n = 5000 and n = 10000. 

n=5000 n=10000 











Mean 


Std 


Mean 


Std 


S Catalan's trisectrix 








20.6684 


0.3691 


20.6966 


0.3051 


L{S) = 20.7846 


L{Sn) 


r 


= 0.5 


20.8722 


0.0697 


20.8289 


0.0479 






r 


= 2 


20.6383 


0.0637 


20.6821 


0.0454 






r 


= 5 


20.6131 


0.0633 


20.6661 


0.0453 


S astroid 


KiS) 






5.5554 


0.1324 


5.6393 


0.1061 


L{S) = 6 


L{Sn) 


r 


= 0.25 


5.7024 


0.0695 


5.7709 


0.0539 






r 


= 0.5 


5.6848 


0.0705 


5.7584 


0.0542 






r 


= 1 


5.6730 


0.0698 


5.7500 


0.0537 
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